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ABSTRACT 

Progress  of  pattern  recognition  as  applied  to 
various  geophysical  signals  is  reviewed.  They 
include  the  teleseismic  signals  for  nuclear 
detection  and  monitoring,  exploration  seismic 
signals  for  oil  prospecting,  marine  seismic 
signals  for  study  of  subsurface  structure  of 
ocean  floor,  seismic  signals  for  earthquake  pre¬ 
diction,  seismic  signals  for  intrusion-detection, 
and  geomagnetic  signals  for  study  of  irregular¬ 
ities  of  the  earth’s  magnetic  fields.  With  the 
exception  of  teloseismic  signal,  research  in 
geophysical  signal  recognition  is  still  in  its 
infancy.  We  expect  to  see  considerable  research 
progress  in  the  next  decade. 

1.  INTRODUCTION 

In  recent  years,  computers  have  played  an  in¬ 
creasingly  important  role  in  gtophysics  especial¬ 
ly  in  seismic  studies  such  as  the  petroleum 
exploration,  mineral  extraction,  nuclear  detec¬ 
tion,  earthquake  research,  marine  seismic  pro¬ 
filing,  and  intrusion  detection,  etc.  Computers 
are  needed  to  process  large  volumes  of  seismic 
data  from  which  useful  information  must  be  ex¬ 
tracted  accurately.  In  addition  to  the  seismic 
data,  geomagnetic  signals  are  also  processed  by 
computer  in  large  quantities.  A  primary  function 
of  computers  has  been  in  filtering  and  storage  of 
the  geophysical  data  while  the  recognition  and 
interpretation  of  the  filtered  results  is  performed 
mainly  by  geophysicists  or  human  operators.  How¬ 
ever  efforts  have  been  made  by  many  researchers  to 
use  pattern  recognition  techniques  in  geophysical 
signal  processing  and  recognition  with  some  promis¬ 
ing  results.  In  this  paper  we  shall  examine  the 
progress  in  geophysical  signal  recognition  and  the 
future  directions  of  such  research.  An  extensive 
list  of  references  provided  as  bibliography  at  the 
end  of  the  paper  clearly  indicates  the  amount  of 
research  activities  in  this  area  in  the  past  ten 
years. 

Pattern  recognition  techniques  considered  in  this 
paper  include  preprocessing,  feature  or  primitive 
selection,  syntactive  analysis,  decision  making 
and  learning  processes.  Among  all  geophysical 
problems,  seismic  discr iminac Ion  clearly  can  be 
formulated  as  a  pattern  recognition  problem  and 
has  been  fairly  extensively  studies  in  this  context. 
For  the  exploration  seismic  data  from  which 
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detection  is  made  of  petroleum,  mineral  or  natural 
gas,  pattern  recognition  can  be  used  as  computer- 
aided  interpretation.  Similarly  for  earthquake 
study  and  seismic  intrusion  detection,  pattern 
recognition  provides  an  additional  dimension  to 
problem  solving.  In  other  geophysical  signals, 
the  need  for  pattern  recognition  is  less  evident. 

As  different  geophysical  signals  require  different 
processing  and  recognition  approaches,  the  recog¬ 
nition  problems  of  various  geophysical  signals 
will  be  considered  individually. 

2.  RECOGNITION  OF  TELESEISMl C  SIGNALS 

Teleseismic  signals  are  the  seismic  waveforms 
resulting  from  nuclear  explosions  and  earthquakes 
which  may  originate  from  unknown  location*. 

Shallow  earthquakes  are  easily  confused  with 
nuclear  explosions.  Seismic  networks  throughout 
the  world  coupled  with  a  priori  knowledge  about 
explosion  and  earthquake  events  have  made  it 
possible  now  to  detect  the  nuclear  explosion  with 
good  accuracy.  Automatic  pattern  recognition  is 
useful  to  aid  or  verify  the  human  interpret  at  ion 
at  the  present  stage  of  development.  Currently 
available  recognition  results  have  not  been  close 
to  error-free  to  make  the  automatic  recognition 
reliable. 

Mathematical  feature  extraction  by  using  dynamic 
spectral  ratio,  complexity,  multivariate  auto¬ 
regressive  parameters  in  conjunction  with  the 
Karhunen-Loeve  expansions  have  been  considered  for 
seismic  discrimination  (e.g.  [ 6 ] ) [i9j [2u] [31j) . 
Features  with  real  geophysical  meaning  appear  to 
be  quite  important  for  seismic  discrimination. 

Examples  are  [19]: 

1.  Earthquakes  produce  approximately  equal  amounts 
of  P  and  S  waves,  while  explosion*  produce  more 

P  waves. 

2.  Earthquakes  have  relatively  ev-cp  foci,  while 
explosions  have  only  shallow  foci. 

3.  Earthquakes  give  anaseismic  ..nd  kataseumic 
first  onsets  while  explosions  give  anaseismic 
first  onset  everywhere. 

4.  The  duration  of  wave  trains  is  shorter  for 
explosions  than  for  earthquakes. 

Unfortunately  the  direct  applic.it ion  of  these 
discriainat ing  features  is  not  possible.  Although 
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explosions  provide  more  P  waves  than  S  waves,  not 
all  earthquakes  follow  the  average  pattern. 

Further  it  is  very  difficult  to  identify  particular 
phases  for  lower  magnitude  earthquakes.  Focal 
depths  hold  a  great  deal  of  promise  in  discrimi¬ 
nation.  However  there  is  still  too  much  error 
associated  with  measurement  to  warrant  its 
exclusive  use.  The  difficulty  in  determining  the 
direction  of  onset,  due  to  noise,  obscures  the 
anaselsmic-kataseisraic  discrimination  feature. 

Again,  the  signal-to-noise  ratio  is  too  small,  for 
small-amplitude  events,  to  discriminate  using  the 
duration  of  the  wave  trains.  An  interesting  dis¬ 
cussion  of  geophysical  features  for  nuclear  detec¬ 
tion  is  given  by  Dolt  [3].  For  automatic  classi¬ 
fication,  the  features  must  be  computable  to  a  good 
accuracy.  This  requires  both  digital  signal 
processing  and  statistical  techniques. 

Figure  1  shows  some  typical  seismic  records  which 
we  have  studied  along  with  their  Fourier  amplitude 
spectra.  The  sampling  rate  is  10  samples/sec. 
with  1200  samples  per  record.  Based  on  our  ex¬ 
tensive  comparative  study  of  various  feature 
selection  criteria  the  spectral  features  when 
accurately  computed  can  perform  the  best  with 
simple  classifier.  For  example,  the  high-reso¬ 
lution  maximum  entropy  spectral  analysis  can  pro¬ 
vide  more  effective  features  than  those  available 
from  the  conventional  modified  periodogram  method 
using  fast  Fourier  transform.  Figure  2  shows  the 
maximum  entropy  spectra  of  two  explosion  and  two 
earthquake  events  for  a  duration  of  200  samples 
from  sample  number  701  to  900.  For  the  explosion 
the  spectral  ratio  defined  as  the  ratio  of  signal 
energy  above  0.5  Hz  to  that  below  0.5  Hz  tends  to 
be  greater  than  1  while  for  the  earthquake  such 
ratio  tends  to  be  less  than  1.  The  nearest 
neighbor  decision  rule  is  used  for  classification. 
For  each  pattern  class,  40  seismic  records  are 
available  for  learning.  A  total  of  243  seismic 
records  which  do  not  include  learning  records  is 
for  classification.  Spectral  ratio  for  200  samples 
from  sample  number  550  to  749  and  spectral  ratio 
for  2C0  samples  from  sample  number  750  to  949  are 
employed  as  two  features  in  classification.  The 
final  classification  result  is  14  errors  out  of 
243  for  94.24%  correct  recobnitlon.  This  is  the 
best  available  result  with  two  features. 

Another  important  feature  extraction  technique  is 
the  short-term  spectral  features.  A  64-point 
window  slides  through  th  1024-point  seismic  record 
generating  64  time  segments  with  48-point  overlap 
between  the  adjacent  segments.  For  each  64-point 
segment,  compute  the  Fourier  transform.  The  first 
32  points  of  the  Fourier  transform  correspond  to 
the  5-Hz  frequency  interval  of  the  seismic  wave. 
Define  -  signal  power  covering  the  first  to 
fifth  spectral  points  (0-0.78  Hz  band),  Pj  "  signal 
power  covering  the  sixth  to  14th  spectral  points 
(0.78-2.18  Hz  band),  and  P3  »  signal  power  covering 
the  15th-32nd  spectral  points  (2.18-5.0  Hz  band). 
Each  F^,  i  *  1,  2,  3  can  be  plotted  as  a  function 
of  time,  i.e.,  the  segment  number.  For  each  curve 
the  ratio  of  the  signal  energies  between  the  first 
32  segments  and  the  second  32  segments  is  called  a 
short-time  spectral  feature.  For  each  seismic 
record,  there  will  thus  be  three  feature  values 


which  form  a  vector.  A  set  of  50  learning 
(training)  records  is  taken  from  each  class.  Using 
the  well  known  Euclidean  distance  nearest-neighbor 
decision  rule  the  best  performance  is  around  80% 
correct  recognition  on  test  set.  For  each  class, 
a  covariance  matrix  Vi(  i  -  1,  2  can  be  determined 
from  the  training  set.  A  modified  procedure  is  to 
compute  the  squared  distance 

(x  -  P(i))  Vt_1(x  -  P(1))  «  d*  ,  i  -  1,2  (1) 

and  choose  the  class  which  provides  the  minimum 
d?.  Here  pCi)  is  a  nearest  neighbor  belonging  to 
the  ith  class  and  x  is  the  three-dimensional 
feature  vector  representing  a  seismic  record.  By 
using  such  weighted  Euclidean  distance  nearest 
neighbor  decision  rule,  a  94.17%  correct  recog¬ 
nition  (13  errors  out  of  223)  is  obtained.  If  the 
Bayes  decision  rule  is  used  under  the  assumption 
of  multivariate  Caussian  density  for  each  class, 
the  performance  is  87.92%  correct  recognition  (27 
errors  out  of  223) .  This  experimental  result 
indicates  that  the  stationary  probability  assump¬ 
tion  of  the  seismic  record  is  not  good.  Other  work 
on  autoregressive  feature  extraction  also  based  on 
the  assumption  of  stationary  random  process  for 
the  seismic  record  has  provided  less  recognition 
performance.  The  short-term  spectral  features 
described  above  have  taken  the  nonstationarity 
into  consideration. 

Among  other  computer  studies,  the  features  based  on 
the  complex  cepstraJ  point  followed  by  Fisher’s 
linear  discriminant  gave  a  94.4%  correct  recog¬ 
nition  results  For  36  seismic  records.  The  auto¬ 
covariance  function  also  can  provide  around  89% 
correct  recognition.  The  results  are  somewhat 
depending  on  the  pre-filtering  operations  used. 

For  a  single  discriminant,  the  ratio  of  surface 
wave  magnitude  M,.  to  the  bodywave  magnitude  M^  is 
most  effective  for  events  with  >  4  as  the 
shallow  earthquakes  typically  produce  more  surface 
waves  than  do  explosions  of  equivalent  energy. 

The  combined  use  of  statistical  and  structural 
features  [7]  In  conjunction  with  an  optimized 
decision  tree  is  a  very  promising  approach. 
Syntactic  method  of  teieseismic  signal  recognition 
has  also  been  examined.  As  a  highly  reliable  and 
fully  automatic  seismic  recognition  system  is  not 
likely  to  be  available  In  the  near  future,  inter¬ 
active  pattern  recognition  is  more  practically 
feasible  and  some  important  effort  has  been  made  in 
this  direction  [21]. 

3.  EXPLORATION'  AND  MARINE  SEISMIC  SIGNALS 

In  the  exploration  seismic  study  detection  is  made 
of  petroleum  or  natural  gas.  In  the  marine 
seismic  study,  an  accurate  mapping  is  made  of  sub¬ 
surface  geological  structure  of  the  ocean  floor. 
Although  the  purposes  are  quite  different,  the  two 
seismic  signal  studies  are  similar.  As  the  new 
petroleum  deposits  are  much  more  difficult  to  dis¬ 
cover  on  land  now,  exploration  of  petroleum  deposit 
in  the  ocean  floor  Is  increasingly  important  and 
the  demand  for  marine  seismic  study  for  oil  ex¬ 
ploration  is  very  high.  For  both  seismic  signals, 
the  major  efforts  made  in  the  past  have  been  in 


digital  filtering  and  estimation  of  subsurface 
parameters.  Pattern  recognition  problems  can  be 
formulated  for  analysis  and  detection  of  the 
seismic  signals.  In  marine  seismic  profiling  [5], 
for  example,  the  physical  parameters  which  describe 
the  structure  of  the  deep  sea  floors  correspond  to 
feature  sets  in  pattern  recognition.  Examples  of 
pattern  classes  are  various  rocks,  sediments, 
mantles  and  oceanic  layers  of  low,  medium  and  high 
porosity.  In  the  exploration  seismic  study  feature 
extraction  andcecision  theoretic  approaches  have 
been  examined  [15] [16].  An  important  recent 
development  is  the  syntactic  approach  [26]  to 
segment  the  one-dimensional  logging  data,  i.e., 
measurements  of  various  physical  rock  parameters 
from  transducers  inside  the  borehole.  The  trans¬ 
ition  patterns  between  signal  blocks  are  described 
by  a  set  of  grammatical  rules. 

As  compared  with  the  teleseismic  signals,  the 
exploration  and  marine  seismic  signals  are  much 
more  stationary  and  thus  both  statistical  and 
syntactic  pattern  analysis  can  be  quite  effective 
for  properly  formulated  recognition  problems. 

4.  MINERAL  EXTRACTION  AND  EARTHQUAKE 
PREDICTION 

The  use  of  seismics  in  mining  exploration  has  been 
restricted  so  far  to  only  some  trial  measurements. 
If  seismic  methods  are  to  become  popular  in  ore 
exploration,  we  can  be  sure  chat  the  signals  of 
interest  will  be  hidden  far  more  "deeply"  than  is 
the  case  for  the  signals  in  oil  propsecting.  The 
highly  developed  reflection  seismology  that 
reveals  subsurface  structure  information  will  be 
useful  to  mineral  extraction.  Pattern  recognition 
techniques  useful  for  peer oleum  exploration 
should  also  be  suitable  for  mineral  extraction. 

Computer-aided  techniques  have  been  quite  helpful 
in  earthquake  prediction.  However  pattern  recog¬ 
nition  techniques  have  not  been  employed  in  earth¬ 
quake  study,  to  our  knowledge;  but  the  potential 
is  there.  For  example,  automatic  processing 
algorithms  for  microcarthquake  data  have  been 
examined  [1],  based  on  the  syntactic  structural 
analysis  of  waveforms  [27] [28].  There  is  always 
ambiguity  associated  with  measuring  the  first 
arrival  time  from  seismograms  whether  it  is  done 
by  a  seismologist  or  by  a  machine  since  the  seismic 
signals  are  of  unknown  shape  and  are  contaminated 
by  noise.  Such  ambiguity  can  be  reduced  by  com¬ 
bining  the  processes  of  picking  arrivals  and  locat¬ 
ing  the  events  in  an  iterative  fashion.  A  seismo¬ 
gram  can  be  described  as  "noise  followed  by 
seismic  signal",  where  both  "noise"  and  "seismic 
signal"  would  have  grammatical  rules  describing 
their  structure.  Current  picking  algorithms  [2] 
use  this  type  of  structural  information  but  in  a 
rather  ad  hoc  manner. 

In  summary,  application  of  pattern  recognition  to 
earthquake  prediction  is  in  its  infancy.  Obviously 
various  structural  and  geophysical  information  must 
be  combined  to  derive  a  feasible  algorithm  for 
automatic  prediction  and  interpretation. 


5.  INTRUSION-DETECTION  USING  SEISMIC  SENSORS 

One  of  the  automatic  intrusion-detection  techniques 
is  to  use  the  seismic  sensors  that  detect  the 
"footstep  signal"  in  the  presence  of  usually 
strong  correlated  noises  due  to  vehicle  motions  in 
the  immediate  vicinity.  In  this  case  the  real 
signal  is  impulse-like  with  unknown  arrival  time 
and  duration.  It  is  necessary  to  employ  noise 
cancelling  algorithms  such  as  the  adaptive  digital 
filtering  and  the  Kalman  filtering  to  suppress  the 
background  noise  so  that  the  footstep  signal  can 
be  enhanced  [10] [11].  Figures  364  show  two 
typical  sections  of  seismic  data  containing  foot¬ 
step  signals,  along  with  the  results  of  adaptive 
digital  filtering  and  adaptive  Kalman  filtering. 

A  major  portion  of  background  noise  has  been 
removed  after  filtering  even  if  the  noise  is  very 
strong.  Each  section  as  shown  has  960  samples  at 
a  sampling  rate  of  200  samples/sec.  For  this  kind 
of  seismic  signal,  spectral  features  are  not  suit¬ 
able  for  detection.  The  recognition  or  detection 
problem  here  is  to  determine  for  each  section  of 
seismic  data  whether  the  footstep  signal  is 
present  or  not.  The  correlation  detection  or 
matched  filter  detection  is  not  useful  here.  A 
simple  peak  detection  with  properly  selected 
threshold  will,  however,  serve  the  purpose.  There 
does  not  appear  to  have  any  "structure"  in  this 
particular  seismic  signal. 

6.  OTHER  GEOPHYSICAL  SIGNALS 

Ceomagnctic  signals  provide  some  information  on 
geophysical  phenomena,  which  is  not  available  from 
the  seismic  data.  The  spectral  analysis  is  the 
main  computer  study  that  has  beer,  performed  with 
the  geomagnetic  signals.  An  example  of  geomagnetic 
signal  is  the  micropulsations  which  arc  often 
associated  wcih  structural  d isturbar.ee-  (or  storms) 
in  the  magnetosphere.  Detecting  a  particular 
micropulsation  event  can  provide  a  diagnosis  of  che 
properties  of  the  magnetosphere.  As  the  magneto¬ 
sphere  changes  because  of  its  varying  interaction 
with  the  solar  wind,  the  properties  of  the  micro¬ 
pulsation  signatures  should  response.  The  recog¬ 
nition  problem  here  is  somewhat  similar  to  the 
detection  of  machine  malfunction.  As  a  largo 
volume  of  geomagnetic  signals  is  received  doily, 
automatic  processing  and  recognition  will  certainly 
be  an  area  of  future  research  study. 

Other  geophysical  signals  include  sunspot  numbers, 
etc.  As  several  different  sensors  are  frequently 
employed  in  examining  a  particular  geophysical 
phenomenon,  how  to  effectively  utilize  information 
from  all  sources  is  a  challenging  problem  in 
recognition  study. 


CONCLUSION 


As  compared  with  recognition  of  speech  and  bio¬ 
medical  signals,  progress  in  geophysical  signal 
recognition  is  much  slower.  There  ire  many 
promising  areas  of  research  as  described  in  this 
paper.  Suggested  future  approaches  are  the 


following. 
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(1)  Understanding  the  geophysical  phenomenon  is 
essential  as  different  geophysical  signals  may 
require  quite  different  recognition  approaches. 

(2)  Signal  processing  is  an  integral  part  of  signal 
recognition  as  a  number  of  distinct  properties  can 
only  be  derived  after  extensive  signal  processing. 

(3)  Fully  automatic  recognition  may  not  be  a 
realistic  goal.  Some  human  interaction  may  be 
necessary  to  complete  the  recognition  task. 

(4)  Information  from  various  geophysical  sensors 
should  be -integrated  end  utilized  to  arrive  at  the 
best  recognition  result. 
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